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SUMMARY 


A study has been made of the accuracy and reliability of current theoreti- 
cal methods of estimating sonic-boom overpressures for level flight at a con- 
stant supersonic Mach number. Theoretical estimation methods have been outlined 
and a numerical evaluation of sonic-boom theory for use on high-speed electronic 
computing machines has been introduced. Estimates based on this theory have 
been compared with available wind-tunnel and flight measurements. The correla- 
tion has been made with wind-tunnel data which incorporate improved data- 
reduction methods. The dependence of sonic-boom overpressure on configuration 
geometry has also been discussed and design methods of minimizing the problem 
have been explored. The results of this review have reaffirmed the conclusion 
that both volume and lift effects contribute to bow-shock overpressures. The 
results have also shown that existing theory provides reasonably accurate esti- 
mates of nominal ground-track overpressures for steady supersonic flight in a 
standard or near-standard atmosphere. 


INTRODUCTION 


The sonic boom, which only a decade ago was an interesting but little- 
recognized and little -understood physical phenomenon, has now emerged as a major 
concern in the operation of present military airplanes and poses one of the most 
serious operational problems to be encountered in the development of commercial 
supersonic transports. In recent years, intensive research efforts treating all 
phases of the problem have served to provide a basic understanding of this phe- 
nomenon. The theoretical studies of references 1 to 3 have resulted in the 
development of estimation methods which have been generally substantiated in 
correlations with the wind-tunnel data of references 4 to 7 and with the flight 
data of references 8 to 13- Effects of atmospheric nonuniformities and airplane 
acceleration and maneuvers have been treated theoretically and experimentally in 
references 14 to 21. Flight data have also provided some knowledge of the 
response of buildings to an imposed sonic-boom overpressure (refs. 22 and 23 ), 
and the psychological reaction of personnel exposed to various overpressure 
levels has also been explored (refs. 23 and 24). 

The purpose of this report is to outline theoretical estimation methods, 
including a numerical evaluation of sonic-boom theory for use on high-speed 


electronic computing machines , and to illustrate the correlation of this theory 
with wind-tunnel and flight measurements. Primary emphasis is placed on tunnel 
data from various sources which have been selected and presented to illustrate, 
in s umm ary form, the more significant findings. Improved data -reduction methods 
and more precise theoretical estimation procedures than were employed in pre- 
vious work have been used in the present correlations of theory with wind-tunnel 
data. Data -re duct ion procedures used herein provide for an adjustment to com- 
pensate for model vibration and other experimental limitations. Throughout this 
report, theoretical estimates are based on area distributions obtained from 
super sonic -area -rule cutting planes. The dependence of sonic-boom overpressure 
on configuration geometry is examined in some detail, and design methods of 
minimizing the problem are explored. Inasmuch as the discussions in this report 
are based, in large part, on wind-tunnel test results, they are necessarily 
restricted to the steady-state case of constant Mach number and altitude. 


SYMBOLS 


A cross-sectional area of airplane or model determined by 

supersonic -area-rule cutting planes having an angle p with 
respect to horizontal 

A(t) nondimensionalized cross-sectional area A/Z^ at nondimen- 

sionalized station t = x/Z 

Afc cross-sectional area at base of airplane or model 


A-e 

Ae,n>A e ,r 


B 

B(t) 

C L 


nondimensionalized effective cross-sectional area due to a com- 
bination of volume and lift effects, A(t) + B(t) 

nondimensionalized effective cross-sectional area A e at non- 
dimensionalized station t = n At and t = r At, respectively 


equivalent cross-sectional area due to lift, 


0 


F L dx 


nondimensionalized equivalent cross-sectional area due to lift 
B/Z^ at nondimensionalized station t = x/Z 


lift coefficient 


Fl 


F(t) 


lifting force per unit length along longitudinal axis of air- 
plane or model 


effective area distribution function, 


_1_ 

2it 




dt 


2 



h 

K r 

l 

M 

n,r 

P 

P 

Pa 

P g 

Ap 

^Pmax 

q. 

S 

t 

At 

W 

x 

AX 

8X 

a 


airplane flight altitude or perpendicular distance from model 
to measuring probe 

reflection factor 

airplane or model reference length 

Mach number 

integers 

reference pressure for a uniform atmosphere (free-stream static 
pressure for wind-tunnel tests) 

mean reference pressure for a nonuniform atmosphere determined 
by method given in reference 18 

atmospheric pressure at airplane altitude 
atmospheric pressure at ground level 

incremental pressure due to flow field of airplane or model 
maximum value of Ap (at bow shock) 

dynamic pressure 
wing planform area 

nondimensionalized distance measured along longitudinal axis 
from airplane nose or model nose, x/z 

incremental nondimensionalized distance along longitudinal axis 
of airplane or model 

airplane weight 

distance measured along longitudinal axis from airplane nose or 
model nose 

distance from point on pressure signature to point where 

pressure-signature curve crosses zero-pressure reference axis 

change in position of bow shock due to vibration 

angle of attack 


P = 



- 1 


3 



7 


ratio of specific heats (1.4 for air) 


4 


Mach angle, sin~l i 

M 


dummy variable of integration measured in same direction and 
using same units as t 

value of t giving largest positive value of integral 
f F(t) dT 


A prime is used to indicate a first derivative and a double prime is used 
to indicate a second derivative with respect to distance. 


NATURE OF THE PHENOMENON 


The nature of the airplane shock field responsible for the sonic -boom 
phenomenon is illustrated in the schematic diagram of figure 1. At supersonic 
speeds the airplane-generated flow field is concentrated within a bow shock and 
a tail shock fanning out from the airplane. When these shock waves reach the 
ground, they are reflected upward. The shock waves traveling with the airplane 
and passing over the ground produce the noise sensed by the observer. Near the 
airplane the pressure signature is quite complex since it contains shock waves 
from the airplane nose, wing-fuselage juncture, engines, tail surfaces, and 
other airplane components. As the distance from the airplane increases, the 
separate shock waves merge and only a bow and a tail shock remain. The result- 
ant N-shape wave signifies the attainment of the so-called far-field conditions. 
At ground level, the incident and reflected signatures are coincident and an 
amplification of the pressure rise occurs. On a hard level surface, a doubling 
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of the pressures takes place. For other surfaces, this reflection factor may 
he somewhat less than 2.0. In free air, in the absence of any reflecting sur- 
faces, this reflection factor may be assumed to be 1.0. 


THEORETICAL CONSIDERATIONS 


The theoretical studies of references 1 to 3 have provided a means of 
estimating the bow-shock pressure rise. In the following equation obtained 
from reference 3 , the bow-shock overpressure directly under the flight path of 
an airplane in level supersonic flight is related to the geometry of the air- 
plane and the flight conditions: 


KrP 1 / 4 


1-197 

\j7 + 1 



( 1 ) 


The function F(t) in equation (l) depends on the longitudinal distribu- 
tion of cross-sectional area and of lift and is defined as follows: 


F(t) 


A. r dt + J_ r B"(t) 

2k J o \J t - t J q \j T _ t 


( 2 ) 


where A (t) represents the second deriv- 
ative of a distribution along the longitu- 
dinal axis of a nondimensionalized airplane 
cross-sectional area determined by 
supersonic -area-rule cutting planes (as 
shown in fig. 2), and B"(t) represents 
the second derivative of a distribution of 
nondimensionalized equivalent area due to 
lift evaluated through an integration of 
the lifting force per unit length along the 
airplane longitudinal axis. Typical air- 
plane distributions are shown in figure 2. 
Since only the pressure field directly 
under the flight path of the airplane is 
being considered in this study, only one 
set of cutting planes having an angle n 
with respect to the horizontal is used. 
Improved accuracy results when the area 
distributions include the increases in 
cross-sectional area due to boundary-layer 
thickness and engine -exhaust effects. 



Figure 2. - Typical distributions used in sonic-boom theory. 
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Inasmuch as the sum of the derivatives is equal to the derivative of the 
sum, equation (2) may he written as 


F(t) 



(3) 


tt 

where A e is the second derivative of an effective cross-sectional area com- 
bining actual cross-sectional area with the equivalent cross-sectional area due 
to lift. This effective cross-sectional area is 


A e = A(t) + B(t) (4) 

The concept of effective cross-sectional area has been used in simplifying the 
numerical method of evaluating sonic -boom characteristics (presented in appen- 
dix A) and has been useful in defining the lower bound of sonic-boom overpres- 
sures discussed in the section entitled "Sonic-Boom Minimization Concepts." 

The length of the positive portion of the pressure signature can be 
expressed by the following equation obtained from reference 3: 



1.19 


7 + 1 M 2 
\[7Tl p 3/4 



(5) 


The slope of the linear portion of the signature may thus be written in the 
* following form, which shows its independence of airplane geometry: 


AP _ k P JL 7 

AX r h M 2 7 + l 


( 6 ) 


The theoretical estimation methods just described have been employed 
throughout the present report. In some of the earlier literature dealing with 
tunnel measurements (refs. 4, 5 > and 7)> these methods were not strictly 
applied. For example, in certain instances a parabolic -body area distribution 
was substituted for the actual distribution and in other instances a distribu- 
tion of normal cross-sectional areas was used instead of that derived from 
supersonic -area-rule cutting planes. 


MODELS, APPARATUS, AND TESTS 


The experimental models and apparatus used in wind-tunnel investigations 
of the sonic-boom phenomenon are unique in several respects. The small size of 
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the models required for an approach to far-field conditions of the flow field 
is perhaps the most unusual feature. For tests conducted in the Langley 4- by 
4-foot supersonic pressure tunnel, model lengths ranging from l/2 to 2 inches 
are required. These small-scale models must he built to extremely small toler- 
ances and must incorporate all the major airplane components. A very sensitive 
pressure -measuring system is also a necessity, since changes in pressure as 
small as 0.02 pound per square foot must be detected in the model flow field. 

An added complication in the design of pressure instrumentation results from 
the relatively large deviations of the reference or free-stream static pressures 
from the nominal value due to time lag in the tunnel control systems. This ref- 
erence pressure may vary as much as 4 pounds per square foot from the nominal 
value in a period of several minutes. Since steady-state tunnel-flow nonuniform- 
ities dictate that the measuring probes or orifices be fixed with respect to the 
tunnel in order to avoid extraneous pressure variations, a means of varying the 
model position during the test must be provided. Because of the vibration of 
the models and apparatus, as well as the presence of a boundary layer and the 
lack of attainment of true far-field conditions, the sharp pressure peaks pre- 
dicted by theory and displayed in signatures measured in flight are not obtained. 
Experimental apparatus and techniques designed to overcome or compensate for 
these difficulties are discussed in the remainder of this section and in appen- 
dix B. 

Two tunnel -apparatus arrangements found to yield satisfactory results in 
sonic-boom tunnel investigations at the Langley Research Center are shown in 
figure 3» In one system, the measurements are made at orifices in a reflection 
plate or boundary -layer bypass plate alined with the tunnel free stream; in the 
other system, the measurements are made by using static -pres sure probes. Both 
systems employ a remote -control sting support for longitudinal positioning of 
the model during the investigations. The apparatus shown in figure 3(b) has 
been used for the more recent tests because it minimizes the boundary-layer 
effect on the measured pressure signatures. 




(a) Plate measurements. (b) Probe measurements. 

Figure 3. - Sketches of typical tunnel apparatus. 
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The schematic diagram of the pressure instrumentation shown in figure 4 is 
applicable for either of the arrangements shown in figure 3- Pressure gages 
having full-scale ranges as low as ±0.05 pound per square inch have been 

employed. The gage measures a pres- 
sure difference between the measuring 
Reference probes or orifices and the reference 

probes or orifices. In order to avoid 
damage to the sensitive gage during 
tunnel startup and shutdown, a safety 
bypass valve has been provided. In 
the design of the instrumentation, 
deviations of the tunnel reference 
pressure from the nominal value have 
to be considered. Since it is only 
a pressure difference that is of con- 
cern, negligible errors result if 
changes in the measured reference pres- 
sure are immediately reflected in 
changes in the measured flow-field 
pressure. In effect, this situation 
Figured- Schematic diagram of pressure instrumentation. can be brought about by equalizing the 

time lag in the instrumentation system 

on the two sides of the gage. Lengths and diameters of tubing are carefully bal- 
anced, and any remaining differences in time lag are compensated through use of 
a small variable -volume device (a bellows). Balancing of the system may be per- 
formed during tunnel evacuation or pressurization prior to a tunnel run. 



Measured 

Adjusted 


A typical measured wind-tunnel pressure signature is shown in figure 5- 
According to theory, supported by flight-test evidence, the pressure signature 
on the ground from an airplane in supersonic flight will have (except at very 

low altitudes) a sharp-peaked N-shape similar to 
Measured that shown by the dashed line. Departures of the 

Ap — . Adjusted measured wind-tunnel pressure signature from a 

] \ true sharp-peaked N-shape wave are caused by the 

presence of near-field effects (double peak) and 
f. \ effects of vibration and probe boundary layer 

Ap i (rounded peaks). In this report the tunnel data 

\ j have been adjusted to compensate for these limita- 

\ f tions simply by extending the linear portion of 

** t the measured curve and forming a right triangle 
ax whose area is equal to the area under the measured 

wave. This adjustment, however, must be applied 
Figure 5.- Typical measured pressure with caution, since it has not been shown to be 

signature. valid where near-field effects predominate and a 

reasonable approach to a far-field signature is 
not evidenced. The maximum value of the measured wind-tunnel bow-shock pres- 
sure rise Ap ma y , as used throughout this study, has been obtained from the 


Figure 5. - Typical measured pressure 
signature. 


adjusted signature. The adjustment (discussed in appendix B) was not applied 
to the data of references 4 and 5; however, a somewhat similar adjustment pro- 
cedure which differed in details of application was used in references 6 and 7- 
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MEASUREMENTS OF THICKNESS -INDUCED PRESSURES 


Wind-tunnel measurements of the flow field surrounding several geometrically 
simple todies (ref. 4) have allowed a study of thickness-induced sonic-boom char- 
acteristics. Representative data 
from that investigation, but now 
incorporating the adjustment dis- 
cussed in appendix B, are presented 
in figures 6 to 8. 


Bow-shock pressure rise 
obtained from adjusted pressure 
signatures for a nonlifting para- 
bolic body of revolution (fineness 
ratio 5) is shown in figure 6 and 
is compared with theory. The theo- 
retical curve presented differs to . 
a small extent from that used in 
reference 4, since in the study 
presented herein the area distribu- 
tion was determined by cutting 
planes inclined at the Mach angle 
rather than by planes normal to 
the body axis. The adjusted data 
display a remarkable agreement with 
the theoretical results both in 
magnitude and rate of decay with 
the distance from the model to the 
measuring probe, even for distances 
as close as 1 body length. Data 
for the parabolic body of rev- 
olution of figure 6 are shown in 
figure 7; in which a pressure 
parameter derived from theory (the 
left-hand side of eq. (l)) is 
employed. A constant value of the 
pressure parameter indicates a 
decay of overpressure with distance 
in accordance with the three- 
quarter-power rule of theory. The 
reflection factor for the boundary- 
layer bypass plate used in these 
tests was assumed to be 2.0. 


AP„ 



Figure 6. - Bow-shock pressure rise for a parabolic body of. 
revolution as a function of distance. M * 2.0; a = 0°. 


Experiment (ref. 4, adj. ) 
Theory 


.16 


.12 


Ap, 


max/h 


3/4 


K r P 


1/4 


.08 


.04 


Figure 7. - Bow-shock pressure-rise parameter for the parabolic 
body of figure 6. 


In figure 8, data are shown for a nonlifting body without axial symmetry. 
Although the overpressure measurements below the model and to the side are sub- 
stantially different at a distance of 1 body length, they become more nearly 
equal at 8 body lengths - a result which indicates an approach to axial symmetry 
of the flow field. Theoretical estimates of pressures to the side of the model 
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Experiment (ref. 4^adj. ) 
Theory 


are identical to those below the model rolled 90° when the supersonic -area-rule 
cutting planes remain at an angle jj. with respect to the horizontal. Area dis- 
tributions used in the theory are shown in the inset sketches. Although the 

shapes of the curves are 
different for the two model 
orientations, the maximum 
areas are identical and the 
theoretical pressure rise 
below the model is only 
slightly greater than that 
to the side. Since these 
differences are extremely 
small, the theoretical 
results presented in this 
figure are represented by a 
single line. 




Figure 8. ■ 


Ap /. N 3/4 
K max/h\ 

P \l 




1/4 


.16 


.12 


.08 


.04 


Bow-shock pressure-rise parameter for a nonaxial 
symmetric body at zero lift. M = 2.0. 


Experiment (ref. 4, adj. ) 
Theory 


its Below 


To the 
side 



Data for a nonlifting 
body that departs radically 
from axial symmetry are 
shown in figure 9- (The 
model is actually a rectan- 
gular wing of aspect ratio 
0.5 with a 12.5-percent- 
thick parabolic -arc section.) 
An examination of these data 
shows that even at 8 body 
lengths there are large dif- 
ferences in the measured 
pressures below and to the 
side of the model and that 
there appears to be little 
tendency toward a further 
approach to axial symmetry 
of the flow field. This 
result is in reasonable 
agreement with the theory, 
which also shows a large 
difference between overpres- 
sures below and to the side 
of the model. These differ- 
ences are directly related 
to the area distributions formed by supersonic -area-rule cutting planes and 
shown on the right-hand side of the figure. The reduced overpressures to the 
side of the model are caused primarily by the increased length and reduced maxi- 
mum cross-sectional area of the corresponding area distribution. These and 
other experimental investigations, both tunnel and flight tests, support the 
validity of the supersonic -area-rule equivalent -body concept in treating 
thickness-induced far-field pressures. 


Figure 9. 


Bow-shock pressure-rise parameter for a rectangular wing at zero lift. 
M-2.0. 
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MEASUREMENTS OF LIFT-INDUCED PRESSURES 


The influence of lift on sonic -boom overpressures may be studied by using 
the data for delta-wing configurations shown in figures 10 and 11. The measure- 
ments of lift-induced pressures shown are taken from the investigation of refer- 
ence 5 but they have now been subjected to the adjustment previously mentioned. 

Bow-shock pressure rise below a 6o° delta wing at angles of attack of 0°, 
5°, and 10° and at M = 2.0 is shown in figure 10. Large differences in over- 
pressure level due to angle of attack (or lift) may be observed in both the 
experimental and theoretical data. These differences may be attributed to dif- 
ferences in the area distributions shown on the right-hand side of the figure. 
The area distributions now include equivalent cross-sectional area due to lift. 


In figure 11, overpressure characteristics for a series of related wing-body 
configurations are presented in a parametric form derived from theoretical rela- 
tionships. Since the wing-body configurations were designed to have the same 
distribution of cross-sectional area and since the delta wings have the same 
equivalent cross-sectional area due to lift at a given value of the parameter 

P, S 
2 L l2 , 


a single curve should serve as the theoretical estimate for the four 


configurations. However, because the cross-sectional area due to the displace- 
ment effect of the assumed laminar boundary layer varies depending upon the model 
wetted area, the theory is presented as a shaded band. Except in the immediate 

vicinity of — Ct = 0, experimental results are in good agreement with theory, 
2 7 2 ' 


both in magnitude and trend. The experimental data shown in this figure are for 
a distance of J2 body lengths; however, the data should apply for greater dis- 
tances since far-field conditions have nearly been achieved, as indicated in fig- 
ure 10 by the small or nonexistent slope of the pressure parameter with distance 
at 32 body lengths. 



Figure 10. - Bow-shock pressure-rise parameter for a lifting delta-wing 
configuration. M ■ 2.0. 


Figure 1L - Bow-shock pressure-rise parameter for a series of 
wing-body configurations as a function of lift parameter. 
M-2.0; h/J-32. 
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CORRELATION OF WIND-TUNNEL AND FLIGHT DATA 


A comparison of tunnel -measured bow-shock pressure rise at M = 2.0 and 
h/l = 50 for a complete supersonic bomber airplane configuration (ref. 6) with 
a theoretical estimate of the pressure rise determined in the manner described 
in appendix A is shown in figure 12. The area distribution used in the theoret- 
ical estimate included the cross-sectional area due to the displacement thickness 
of a laminar boundary layer. The figure shows close agreement between wind- 
tunnel measurements and theoretical overpressure values. 


Flight-test results for the supersonic bomber airplane are reported in ref- 
erence 12. These tests have produced the most extensive and most self-consistent 
data yet recorded for ground measurements of the sonic boom created by an air- 
plane. The instrumentation employed is noteworthy in that it faithfully repro- 
duced the entire pressure signature and not just the peaks as in previous tests 
(for example, refs. 9 to ll) . The data are of particular interest since the 
altitude range extends to 75 >000 feet where lift effects are important. Corre- 
lation of these flight data with theoretical estimates in a manner similar to 
that used for the correlation of the tunnel data is possible for this same air- 
plane. However, for use in equation (l), it is first necessary to define a ref- 
erence pressure which accounts for the variation of atmospheric pressure and 
temperature between the airplane and the ground. In most previous work, a crude 
approximation, the geometric mean of the atmospheric pressure at altitude and 


that on the ground 




has been used. A thorough study of the effects of 


the atmosphere given in reference 18 provides an evaluation of shock strength as 
the shock propagates through the layers of an assumed stratified atmosphere. 

From such information for a standard atmosphere supplied by the authors of ref- 
erence 18 , a mean or effective reference pressure p has been evaluated and 
the ratio P /l/P a P e i s presented in figure 13 . Note that at* the altitudes nor- 
mally associated with supersonic 
flight, p is substantially 
greater than JP a Pg* 




Experiment (ret. 6, adj. ) 
Theory 



Figure 12. - Bow-shock pressure-rise parameter for a wind-tunnel model of a 
supersonic bomber airplane. M - 2.0; h/l » 50. 


Figure 13. - Ratio of mean reference pressure 
(determined by method of ref. 18) to geometric 
mean. 
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With the use of the mean reference pressure p, the flight data for the 
supersonic bomber airplane (ref. 12) may be reduced to the parametric form shown 
in figure 14. Each data point repre- 
sents the average pressure rise 
recorded by several microphones in 
use during a given overpass. These 
microphones were spaced along the 
airplane ground track over a distance 
of about 4 miles. A reflection fac- 
tor of 2.0 was used in evaluating the 
pres sure -rise parameter since meas- 
urements during the test indicated 
that the dry lakebed over which the 
flights were conducted acted as a 
nearly perfect reflection plane. The 
flight data show overpressures 
slightly higher than the theoretical 
estimates in which an increase in 
cross-sectional area due to a turbu- 
lent boundary layer has been included. 

On the average , the measured results 
fall about 5 percent higher than the 
theoretical estimates; however, the 
variation with lift parameter is in 
close agreement with the theory. The 
tunnel data of figure- 12 
and the flight data of fig- 
ure 14 provide strong evi- 
dence that lift effects are 
well predicted with 
present-day theoretical 
estimation methods. In 
addition, these data tend 
to indicate that atmos- 
pheric effects may be bet- 
ter accounted for by using 
the methods of reference l8 
than by using previous 
methods . 



Figure 14. - Flight measurements of bow-shock pressure-rise 
parameter for a supersonic bomber airplane. M - 1.5 to 2.0. 


Theory 

Flight (ref. 12) 



A direct comparison of 
pressure signatures meas- 
ured in the flight test 
with a signature measured 
in the tunnel test is shown 
in figure 15- This compar- 
ison has been made possible 
through a correspondence in 
Mach number and lift coef- 
ficient for a flight data point and a tunnel data point. Representative signa- 
tures from four of the microphones in use during the overpass (including the max- 
imum and minimum pressure peaks) are compared with the tunnel measurement and 


Figure 15. - Comparison of flight and tunnel measurements of pressure signature 
of a supersonic bomber airplane. M - 2.0; | C^-| ■ 0.0135. 
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pq 

6 \ 


Ap 

K max 

Ib/sq ft 


Flight (ref. 12) 
Tunnel (ref. 6, adj.) 
Theory 


x. 


with theory. The dashed line, 
which represents an adjusted 
tunnel signature as defined in 
appendix B, displays a reason- 
able agreement with the flight 
signatures . 

A more familiar form for 
the presentation of flight 
sonic -boom data than the para- 
metric form used in figure l4 is 
given in figure 1 6. Average 
values of the maximum ground- 
measured overpressures are plot- 
ted as a function of altitude. 

The estimates, based on theory 
and on extrapolated tunnel data, 
are shown as a band in order to 
account for variations in flight 

Mach number and airplane weight at a given altitude. The range of Mach numbers 
from 1.5 to 2.0 and the range of airplane weights from 62,000 to 92,000 pounds 
account for a much narrower band than is shown by the scatter of the data. 
Nevertheless, the agreement of the flight data with the estimates is generally 
good and indicates, to some extent, the degree of confidence which may be placed 
in further estimates of nominal ground overpressure using these methods. 


30 


40 


50 60 

Altitude, ft 


70 


lx l(f 


Figure 16. - Flight measurements of ground overpressures for a 
supersonic bomber airplane. M - 1.5 to 2.0. 


SONIC -BOOM MINIMIZATION CONCEPTS 



The theory of references 1 to 5 has made possible the definition of a lower 
bound of sonic-boom overpressure, which has been discussed in references 25 and 
26. As shown previously, sonic-boom strength depends on an effective area dis- 
tribution in which both volume and lift components are combined. An example of 
an effective-area-distribution curve is shown in figure 17. Note that the value 
of A e at the base of the airplane is fixed by the 
airplane base area (including boundary-layer and 
engine-exhaust areas) and by the flight conditions 
of Mach number and lift coefficient. Although the 
sonic -boom strength parameter has been found to 
depend primarily on the value of the effective 
cross-sectional area at the base, it also depends on 
the shape of the complete A e curve. In refer- 
ence 25 the shape of the area -distribution curve 
yielding a minimum sonic boom was shown to be given 
by a function in which the area is proportional to 
the square root of the distance from the nose except 
in the immediate vicinity of the airplane nose. 

Such a curve is shown by the dashed line of fig- 
ure 17. As shown in reference 26, a lower bound of 
attainable sonic-boom overpressure that depends only 



Typical 

Lower bound 


3 ' r 1 
2 C L ,2 


Figure 17. - Lower-bound effective- 
area distribution. 


on the airplane length, weight, base area, and flight conditions can he written 
in simplified form as 


4Pmax/ hf /4 

P \l) 

Krp 1 / 4 



(7) 


Although this form of the equation follows that used in reference 2 6 , the shape 
of the optimum-area curve and the resultant shape factor of 0 - 5 4 were obtained 
with the use of reference 25- The lower-bound expression neglects any minimum- 
volume restrictions and thus is inapplicable near zero lift coefficient, as 
explained in reference 26. This omission has not proved to be serious in the 
studies made to date. When volume restrictions become necessary, they may be 
included, as was done in reference 26. In the limiting case when A -5 is zero, 
the lower-bound expression may be written as 


A Pmax ^ 3/4 

K r pV 4 


°-5M 5 C L To = 0.5M & 


3 w 


2 q 


(8) 


Experiment 
(ref. 7, adj. ) 
a 


Theory 




The lower-bound concept not only sets limits on attainable overpressures 
but also suggests design methods of approaching these limiting values. Theoret- 
ically, for a selected flight condition (a design point), it should be possible 
to design a configuration to 
approach the sonic -boom minimization 
requirements . Some experimental 
data (ref. 7 ) believed to be appli- 
cable in connection with these con- 
cepts are shown in figure 18. Meas- 
ured and theoretical overpressures 
in parametric form have been plotted 
against the lift parameter for two 
wing-body models. The model with 
the wing in the rearward location 
theoretically approaches the lower 
bound even though it was not 
designed strictly in accordance with 
the concepts previously discussed. 

Effective-area-distribution curves 
for an assumed design point 


■'\r 


^2 Cl 2 = 0.008 6j are shown on the 

right-hand side of the figure. 

These curves illustrate graphically 
that the area distribution of the 
rear-wing model more closely 




Figure 18. - Comparison of sonic-boom characteristics of 
configurations approaching, to different degrees, the 
requirements for lower bound. M - 2.0; h/I - 50. 
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approaches that of the lover bound. The experimental results confirm the trends 
predicted by theory. Reference 27 provides a study of the effects of the neces- 
sary compromises with airplane drag on the attainment of sonic-boom lower bound. 



Figure 19. - Comparison of sonic-boom characteristics of two 
supersonic transport configurations. M ■ 1.4 and 2.0; 
h II = 50. 


In order to provide an illus- 
tration of configuration effects, 
a comparison of the sonic -boom 
characteristics of two transport 
configuration models is shown in 
figure 19. Both theoretical and 
experimental wind-tunnel data are 
shown and are compared with a lower- 
bound curve for which the base area 
A b is assumed to be zero. Inas- 
much as theoretical differences in 
the sonic-boom characteristics at 
the two Mach numbers are small, 
only a single curve is shown for 
each configuration. Cross-sectional 
areas used in the theory include the 
area within the estimated displace- 
ment thickness of a laminar boundary 
layer. The lower overpressures for 
the arrow-wing design may be attrib- 
uted to the reduced base area and 
to the smoother area- and lift- 


distribution curves. Examples of 

effective-area-distribution curves for an arbitrarily selected value of the lift 
parameter are also shown in the figure. These data illustrate the significant 
effect of configuration geometry on sonic -boom strength. 


CONCLUDING REMARKS 


A review and analysis of wind-tunnel and flight sonic-boom data incorpora- 
ting recent developments in wind-tunnel data-reduction methods has reaffirmed 
the conclusion that both volume and lift effects contribute to bow-shock over- 
pressures. The study has also shown that existing theory provides reasonably 
accurate estimates of nominal ground-track boom overpressures for steady flight 
in a standard or near-standard atmosphere. It has been indicated, however, that 
design considerations based on developments of the theory will be only partly 
effective in minimizing the sonic -boom problem. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 7, 1964. 
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APPENDIX A 


A NUMERICAL EVALUATION OF SONIC -BOOM THEORY FOR USE ON 
HIGH-SPEED ELECTRONIC COMPUTING MACHINES 


Theoretical estimates "based on the work of references 1 to 3 and used in 
this report for correlation with the experimental data have been evaluated by 
using a numerical procedure. Because of the 
nature of the curve-fitting technique employed, 
the area distributions are made to be smooth 
(that is, having no discontinuity in the first 
derivative ) . Thus , rather than by employing 
more rigorous solution suggested in refer- 
ence 2, discontinuities in the first deriv- 
ative are accounted for by assuming the 
changes to occur over finite distances. 



As discussed previously, configuration 
sonic-boom characteristics are directly 
dependent on an effective-area-distribution 
curve A e formed by a direct addition of 
actual area and equivalent area due to lift. 

As shown in figure Al, the A e curve may be 
approximated by a series of parabolic arcs 
having a first derivative composed of contin- 
uous straight-line segments and a second 
derivative composed of a step or pulse func- 
tion. The integral involved in the F(t) 
function can be evaluated quite easily when 
A e is a constant; and by superposition, a 

complete F(t) curve may be built up corre- 
sponding to the Ae pulse distribution. An 
integration of the F(t) function to the 
point t q (cross-hatched area in fig. Al) is 

then used in evaluating the right-hand side of 
equation (l). The degree of approximation of 
the A e curve can be improved by increasing 

the number of pulses used. 


A e 






i 

t 



Figure Al. - Numerical method of determining 
sonic-boom characteristics. 


The effective nondimensionalized cross-sectional area may be expressed as 



dx 


dx 
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or 


A e 


A(t) + | C L 


s 

Mil 


B(l) 


(Al) 


If it is assumed, that local lifting pressures are directly proportional to the 
total lift coefficient, the quantity in brackets in equation (Al) is a constant 
for any airplane station. The inputs to the program are thus tabulations of 
A(t) and B(t)/B(l) as functions of equally spaced nondimensionalized airplane 
stations. Equation (Al) then allows an evaluation and tabulation of the effec- 


tive area distribution A e 


for selected values of the lift parameter 


P c 

2 C L 



A series of parabolas, fitted to these points so that the resultant curve has 
no discontinuities in slope, may be expressed as follows: 




A e,l 2 

(At) 2 


(0 < t < At) 


Ae ,2 = A e 


, , ( A e,2 - A e,l) - ( A e,l) At , 

(At f 

+ ( Ae ' 3 - Ae,2) - ( A e,2)< 


At) c 


(At < t < 2At) 


> (A2) 


(2At < t < 3At) 


A e,3 = A e , 2 + Ai j8 (t - 2At) + - 2At) 2 

(At) 2 

, r . , -l (A«jn - A e,n-l) - ( A e,n-l)At r -.2 . \ 

A e,n = A« jn _i + A e ,n-l[t - (n - l)At| + A — — — i L — Jt - (n - l)Atj ((n - l)At < t < nAt) 


(At) 


The second derivative is 


A e,l o A e,l 

(At) 2 


A e,2 = 2 g ( A e,2 - A e,l) ' a e,l 

(At) 2 ^ 


(0 < t < At) 


(At < t < 2At) 


> (A3) 


A 'e,3 = ^l( A e,3 ‘ A e, 2 ) - 2 ( A e,l + A e, 2 ) (2At < t < 3At) 


a" - 2 

A e,n - 


(At) c 


/ \ r=n-l 

( A e,n - ke,n-l) - 2 , T f( n - l)At < t < nAt) 

' ' r=l 
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The F(t) function then becomes 


f(t) = \ Ae^fT 


F(l-) = i 


F(t) = i 


a + (Ae j2 - ^,i)/rrst 

A^,n/T + (a’’, 2 ' Ae A )\|r - At + (a", 5 - 4,2 )|/7T 


2At 


(0 < t < At) 
(At < t < 2At) 
( 2At < t < 3At) 


> (A4) 


F(t) = - Ae,r-1 \/ T - (r - l)At 

« r=l ' 7 


((n - l)At < t < nAt) 


and the integral of the F(t) function may be written as 


•>T=At 


F(r)dr = A" x (At) 


' T=2At 


F(T)dT 


3« 


pT=3At 

/ F(-r)dT = — 

'o 3,r 


T 


^T=nAt 


F(T)dr = 


3« 


jAe,;L(2At) 5 / 2 + (4,2 - Ae,l)(At) 3 / 2 

A^jOAt) 3 / 2 + (a’^ 2 - A e,l)(2At) 3 / 2 + (a’^j - Ag j2 )( At 

3 /2 

^Ae,r ■ A ejr _i)[](n - r + l)At] 


) (A5) 


r=l 


The area -distribution curve-fitting technique employed produces solutions for the 
derivatives, the F(t) function, and the integral of F(t) that oscillate from 
point to point. From comparisons of numerical solutions with certain analytical 
solutions for simple bodies, these oscillations have been found to center on the 
analytical solutions. Improved accuracy results when the integral of F(t) is 
averaged as follows : 
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Finally, the overpressure is found by selecting the maximum value of equa- 
tion (A6) and substituting it into the following equation: 



1.197 


F 1 b 


'T=nAt 


F(T)dr 


max 


(A7) 


Equations (Al) to (A7) can be readily adapted for use on high-speed elec- 
tronic computing machines in a numerical evaluation of sonic-boom theory. The 
computational program may be summarized as follows. The inputs are tabulations 
of the nondimensionalized area A(t) and the integrated lift distribution ratio 
B(t)/B(l) as functions of equally spaced nondimensionalized airplane stations. 
Equation (Al) then allows the determination of a table of effective cross- 


sectional areas for preselected values of the lift parameter 


P r _S_ 

2 L l2 - 


The 


second-derivative step function is then generated by using equation (A3). Tab- 
ulated values of A e are used in evaluating equation (a 6), the maximum value of 


this integral being selected and used in equation (A7) to evaluate the sonic- 
boom characteristics of the configuration. For airplane configurations employing 
camber, the loading distribution at zero lift may be taken into account by using 

a 

a modified area distribution. In this case — — is replaced by 

l 2 

a b r x 

— — + — / F t dx where F L is the lifting force per unit length at zero 

l 2 2ql 2 J 0 


lift. 
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ADJUSTMENTS OF WIND-TUNNEL MEASUREMENTS OF BOW-SHOCK STRENGTH 


TO COMPENSATE FOR EXPERIMENTAL LIMITATIONS 


A number of experimental difficulties arise in attempting to measure 
within the confines of a wind tunnel the pressure signatures of the necessarily 
small models and in attempting to extend the results to apply to full-size air- 
planes at flight altitudes. The necessity of attaining or approaching far-field 
conditions, where the pressure signature assumes a characteristic N-shape, 
requires that tunnel models be extremely small. Even with models as small as 
those employed in the investigations of this report, an approach to far-field 
conditions is not achieved in all cases. It does not appear to be practical to 
reduce further the model size because of construction difficulties and because 
vibrations of models, probes, and support apparatus introduce changes in the 
shape of the pressure signature and in the magnitude of the pressure rise, which 
become progressively more pronounced as model size is decreased. The presence 
of a boundary layer on the measuring probe also introduces changes in the shape 
of the signature and in the magnitude of the pressure rise, which are dependent 


on model size. Another result of 
decreased model size is the increase in 
relative importance of the increment in 
effective cross-sectional area due to 
model boundary layer. 

With a compromise model size, the 
experimental limitations in attaining a 
far-field N-shape pressure signature are 
always present to some degree and are 
occasionally large. Thus, a method of 
interpreting the results and compensating 
for these limitations becomes necessary. 
The following discussion explores these 
problems and describes the method used to 
adjust the wind-tunnel measurements of 
bow-shock strength. 

The failure to achieve a classical 
N-shape wave in tunnel tests is due in 
part to the fact that in many cases the 
pressure signatures are in the transition 
region from near-field to far-field condi 
tions, as shown by the data of figure Bl. 
These measurements made at M = 2.0 for 
the canard transport configuration shown 
in figure 19 have been plotted in a para- 
metric form suggested by theoretical con- 
siderations. According to the theory, 



Figure Bl. - Transition of pressure signature from near-field 
to far-field form. 
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when far-field conditions (an N-wave) have "been achieved, signatures plotted in 
this form remain identical as distance is increased. The near-field shape of 
the pressure signature at h/z = 12.5 is evidenced by the presence of two dis- 
tinct pulses in the region of the bow shock. These apparently are the separate 
shocks from the fuselage nose and from the wing -body juncture. At a distance of 
50 body lengths, the pulses have merged and an N-shape wave is approached. It 
has been noted that even for quite complex signatures, a linear portion of the 
pressure signature develops and the slope closely agrees with that estimated by 
far-field theory. By accepting the premise that, during this transition (as far- 
field conditions are approached) , the impulse area under the bow-shock portion of 
the signature attenuates with distance in a manner similar to that for a fully 
developed N-wave, an attempt may be made to define the pressure signature that 
would exist if far-field conditions were established. The adjusted signature 
may be determined as illustrated in figure B1 simply by extending the linear por- 
tion of the measured signature forward so that a right triangle is formed whose 
area is equal to the area under the measured curve. Because of inexactness in 
the assumptions, the adjustment cannot be rigorously correct; however, a practi- 
cal test would appear to be met when adjusted signatures plotted in the form used 
in figure B1 remain constant as distance is increased. The remaining discrepan- 
cies between the measured signature at h/Z = 50 and the sharp-peaked N-shape 
signature which would be expected in the far field are believed to be caused by 
vibrations of the model and support apparatus as well as by boundary-layer 
effects . 

In order to study the influence of vibration, consider a completely steady 
model in uniform supersonic flow and an ideal pressure-sensing system with a 
probe at a distance large enough to enable a true far-field N-wave to be 
recorded, as represented by the long-dash line in figure B2. Suppose that the 
model (or the measuring probe) undergoes a constant -amplitude vibratory motion 
represented by the inset sketch. The N-wave will then occupy successive posi- 
tions at equal time increments as indicated by the short-dash lines on the pres- 
sure signature plot. At a given longitudinal probe location, a highly damped 
measuring system such as the one used for these tests would register a time 
average of the pressures imposed on it. When a range of probe locations is con- 
sidered, the measured pressure signatures with a constant -amplitude vibrating 
system take on the appearance of the solid-line curve. This curve does not 
resemble the actual wind-tunnel data, but it is not likely that tunnel vibration 
is confined to the single amplitude shown in figure B2. When a varying amplitude 
is considered, the resulting pressure signature assumes the characteristics of 
that shown in figure B3. The assumed amplitude -time relationship is shown in the 
inset sketch. The resulting signature now resembles those obtained from actual 
tunnel measurements . 



Figure B2. - Effect of constant-amplitude vibration. Figure B3. - Effect of varying-amplitude vibration. 
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In both figure B2 and figure B3, note that the areas under the curves are 
almost unchanged from the steady to the vibrating condition. Also note that the 
middle portion of the signature remains unaffected provided the amplitude of the 
vibration is less than the length of the signature. These observations may now 
be utilized in an attempt to adjust the measured data to provide an estimate of 
the pressure signature in the absence of vibration. This adjustment may be 
accomplished by extending the linear portion of the measured signature forward 
so that a right triangle is formed whose area is equal to the area under the 
measured curve. Since this procedure is identical to that previously discussed 
in the compensation for the presence of near-field pres sure -signature character- 
istics, one adjustment will suffice for both deficiencies. 

The foregoing discussion of vibration effects was considered to be independ- 
ent of possible viscous effects. The probe boundary layer, however, is a signif- 
icant factor in the sensing of static -pressure changes across shock waves. The 
imposition of shock-wave pressure gradients on boundary layers of pressure- 
sensing instruments generally produces flow distortions which can be sensed both 
upstream and downstream of shock locations. This condition effectively results 
in tendencies for instrument -sensed pressure changes across shock waves to be 
less abrupt than pressure discontinuities across the shock waves in the absence 
of instruments. Such effects of boundary layer, as well as effects of vibration, 
in spreading and rounding off shock-wave pressure signatures are approximately 
accounted for by the previously described technique for adjusting wind-tunnel 
pressure measurements. The applicability of the adjustment technique may be 
uncertain, however, if the pressure-sensing arrangements are different from those 
employed in these investigations. 
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